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INTRODUCTION 


The  change  in  velocity  of  an  acoustic  wave  in  a  solid  due  to  stress  is 
known  as  the  acoustoelastic  effect.  Hughes  and  Kelly  (ref  1)  as  well  as  Bach 
and  Askegaard  (ref  2)  have  derived  expressions  for  the  velocities  of  plane 
acoustic  waves  in  homogeneously  stressed,  isotropic,  and  homogeneous  solids. 
These  expressions  show,  to  a  first-order  approximation,  that  the  relative  change 
in  velocity  (Av/v)  is  proportional  to  the  uniaxial  stress  or  a  linear  com¬ 
bination  of  the  triaxial  principal  stresses,  with  coefficients  that  are  func¬ 
tions  of  the  second-  and  third-order  elastic  constants.  The  calculations, 
unfortunately,  are  not  readily  extendable  to  Rayleigh  waves  or  inhomogeneous 
stress  situations. 

Hayes  and  Rivlin  (ref  3)  calculated  the  acoustoelastic  effect  for  Rayleigh 
waves  propagating  on  the  surface  of  a  uniformly  stressed  material.  These  calcu¬ 
lations  were  later  extended  by  Hirao,  Fukuoka,  and  Hori  (ref  4)  to  one  par¬ 
ticular  configuration  of  Rayleigh  wave  on  an  inhomogeneous ly  stressed  medium. 
Again,  the  calculations  have  the  disadvantage  that  they  are  not  directly 
applicable  to  situations  involving  other  configurations  of  Rayleigh  wave  and 
inhomogeneous  stress,  including  the  one  considered  in  this  report. 

The  perturbation  theory  for  acoustoelastic  effects  as  recently  developed  by 
Husson  and  Kino  (refs  5,6),  however,  is  quite  general  and  can  be  applied  in  a 
straightforward  way  to  various  configurations  of  acoustic  waves  and  homogeneous 
or  inhomogeneous  stress  states.  We  have  applied  this  theory  to  the  Rayleigh 
wave  situation  as  encountered  in  our  work. 

References  are  listed  at  the  end  of  this  report. 


THEORY 


The  starting  point  for  applying  this  perturbation  theory  to  Rayleigh  waves 
is  the  corrected  form  of  Eq.  (20)  from  Reference  6: 

3  b 

s<t>  =  ~  2p  f  { (2i+X)  [A(a2)  +  B(a2>  +  C(a2>]  +  (X+m)D(a2>  +  mE(a2)} 

<*r  y  oam 

3bi  n 

+  3ai  (x+2m-n)c(a2)  -  2  [D(a2)  +  E(a2)]} 

3b2 

+  ai~  {(2X+6p+4m)A(a2)  +  p[2D(a2)  +  E(a2)]} 

3b3 

+  { (2X+6pt-4m)B(a2)  +  p[2D(a2)  +  E(a2)]})dV  (1) 

(Here  and  throughout  the  report  we  employ  the  Einstein  summation  convention.) 

In  this  equation,  6<t>  is  the  phase  shift  experienced  by  the  Rayleigh  wave. 
dbm/3am  are  the  derivatives  of  the  static  displacements  (static  strains).  X,  u, 
1,  m,  and  n  are  the  Lam6  and  Murnaghan  constants  (second-  and  third-order 
elastic  constants).  The  a,  are  the  coordinate  axes,  with  a3  in  the  direction  of 
propagation  and  a2  normal  to  the  surface  and  directed  inward,  w  is  the  angular 
frequency.  P  is  the  average  power  flow  given  by  (ref  6): 
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(2) 


where  p  is  the  density  of  the  undeformed  material,  vr  is  the  Rayleigh  wave 
velocity,  and  u,  are  the  particle  displacements  for  the  Rayleigh  wave.  A,  B,  C, 
D,  and  E  are  functions  of  a2  as  determined  from  the  following: 


The  particle  displacements  u-j  for  a  Rayleigh  wave  propagating  along  83  are 
given  by  (refs  6,7,8): 

ui  =  0 

2  r@"aSa2  .  #»-aLa2)e_i(/Sa3 -«t) 
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The  shear,  longitudinal,  and  Rayleigh  wave  velocities  (vg,  vL,  and  vR)  are  in 
turn  given  by  (refs  7,8): 


vs  = 


VL 


X  +_2p 
P*' 


,0.87  +  1.12V. 

UR  ■ 


with  v,  Poisson's  ratio,  given  by 
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(7) 


After  evaluating  the  integrands  in  Eqs.  (1)  and  (2)  over  a  unit  length 
along  a^,  a  length  L  along  83,  and  from  0  to  ®  along  83,  we  calculate  the  rela¬ 
tive  change  in  velocity,  Av/v,  using  (ref  6): 


Av  _  9b3  5<f>vR 

v  883  gjL 


(8! 


Of  course  to  perform  the  integrations,  one  needs  to  assume  a  form  for  the  static 
strains  consistent  with  a  load  applied  as  shown  in  Figure  1.  We  assume  that  the 
normal  stresses  vary  linearly  with  depth  a 2 : 

ffl  =  0  *  °2  ~  (H2a2)ao  »  a3  =  (l+H332)ao  (9) 


where  a0  is  the  surface  stress  in  the  83  direction,  and  the  constants  H2  and  H3 
must  be  determined  by  applying  the  principles  of  elastic  theory  to  the  geometry 
of  our  specimen  and  load.  Application  of  Hooke's  law  will  then  yield  the  static 
strains  to  use  in  Eq.  (1). 

The  assumptions  of  Eq.  (9)  lead  to 
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where  the  are  complicated  functions  of  the  second-  and  third-order  elastic 
constants  and  the  density.  They  are  best  evaluated  numerically  rather  than  ana¬ 
lytically.  The  acoustoelastic  constant  for  Rayleigh  waves  in  this  situation 
would  then  be  equal  to  the  term  in  parentheses. 


CALCULATIONS 

For  the  type  of  steel  used  in  this  study  (ASTM  A723/MIL-S-46119A) ,  the 
elastic  constants  and  density  have  been  found  by  Scholz  and  Frankel  (ref  9): 

A  =  110.3  GPa 
jj  =  79.9  GPa 
1  =  -388  GPa 
m  =  -624  GPa 
n  =  -747  GPa 


p  =  7.84  g/cm3 

Using  these  values,  we  obtain  yx  =  -0.00411  1/GPa,  y2  =  -84.7  m/(GPa-s),  and  y3 
=  -79.7  m/(GPa-s) . 

To  determine  H2  and  H3,  we  first  obtain  the  stresses  in  the  ring  as  a 
function  of  radial  position  r  using  elastic  theory  (ref  10): 
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,v 


-F/D  f  -a2b2  r  1 

a2 ' r '  =  1-7a*+b*)r3  '  a2“+“b2  +  rJ 

ln(b/a)  - 

b2  +  a2 


2(a+b)F/D  . (b/a) 2 ln(r/b)  -  ln(r/a)  +  (b/r ) 2 ln(b/a) 

l2abTn(b7aI]2  1  (b7a)2"--I  .  J 


b2  -  a2 


c»3(r)  = 


1n(b/a)  - 


)  .  a2b2  3r  1 

b2~-~a2  L(i2+b2Tr®  "  a2"+“b2  +  rJ 


b2  +  a2 


.  _ ?^a+b]F/D . .  .  ib/al21n{r/b2  -  ln{r/a)  -  {b/r)21n{b/a) 

[2ob1n(b/a) ] 2  u  """  (b7a)2_-'l  J  Ul* 

b2  -  a2  -  - 

b2  -  a2 

F  is  the  force  compressing  the  ring  at  the  split,  D  is  the  thickness  of  the  ring 
in  the  axial  direction,  and  a  and  b  are  the  inner  and  outer  radii,  respectively. 
This  solution  is  correct  only  in  the  region  opposite  the  split.  By  dividing  the 
derivatives  of  o2  and  <73  by  a3  at  r  *  a,  we  obtain  H2  *  12.42  1/m  and  H3  = 

-497.8  1/m  (F  and  0  drop  out). 

Using  the  above  values  of  the  y's  and  H's  in  Eq.  (10),  we  obtain 

Av  Bi 

--  =  <B0  ♦  f~)o0  (12) 

where  B0  =  -0.00411  1/GPa,  Bj  =  0.00615  MHz/GPa,  and  f  is  the  frequency  in  MHz. 
Figure  2  shows  the  theoretical  acoustoelastic  constant  versus  frequency. 

EXPERIMENTAL  METHODS 

The  specimen  was  a  split  ring  of  ASTM  A723  steel  (MIL-S-46119A)  machined  to 
the  dimensions  shown  in  Figure  1.  An  inhomogeneous  stress  was  applied  using  a 
compressive  load  at  the  split.  We  used  three  longitudinal  wave  transducers 
mounted  on  plastic  wedges  cut  to  the  critical  angle  to  couple  a  longitudinal 
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wave  in  the  wedge  to  a  Rayleigh  wave  on  the  surface.  All  three  wedges  were 
bonded  to  the  steel  with  silicone  rubber  (GE  RTV) .  The  Rayleigh  wave  was  intro¬ 
duced  at  the  first  wedge  and  received  at  the  second  and  third  ones. 

Velocity  data  were  taken  using  a  computer-controlled  measurement  system 
(Matec  Instruments,  Model  MBS-8000)  based  on  phase  detection  methods  developed 
by  Peterson  (ref  11).  The  block  diagram  of  Figure  3  shows  the  system's  main 
features.  The  measurement  technique  involves  interactive  automatic  control  of 
i  the  frequency  and  measurement  of  phase  relationships. 

From  the  phase  detectors  shown  in  the  figure,  the  computer  receives  signals 
proportional  to  the  sine  and  cosine  of  each  received  pulse's  phase  with  respect 
to  the  original  reference  wave.  The  computer  can  then  calculate  the  amplitude 
and  phase  of  each  received  pulse.  Of  necessity,  the  phase  is  calculated  only  as 
an  angle  between  -71  and  +tt.  Using  an  algorithm  that  interactively  varies  the 
frequency  slightly  and  measures  the  corresponding  phase  shifts,  the  system 
calculates  the  transit  time  of  the  acoustic  wave  and  the  tvie  difference  between 
receipt  of  the  pulse  at  the  two  receiving  transducers.  The  system  also  calcu¬ 
lates  changes  in  that  time  difference,  due  to  the  application  of  various  loads, 
by  measuring  the  corresponding  phase  shifts. 

For  ease  of  discussion,  let  the  transit  time  from  the  sending  transducer  to 
the  first  receiving  transducer  be  Tj,  and  to  the  second  receiving  transducer, 

)  T2.  We  then  measured  the  time  change  AT  =  AT2  -  ATj  versus  applied  load.  By 

measuring  the  strain  on  the  outer  surface  of  the  ring  (Figure  1)  and  using  the 
principles  of  elastic  theory  (ref  10),  we  calculated  the  applied  strain  on  the 
inner  surface  between  the  two  receiving  transducers.  The  relative  change  in 

velocity  Av/v  can  then  be  calculated  using 

Av  AT 

v*  =  £ID  ’  I” 


where  £iq  is  the  strain  at  the  inner  diameter  and  T0  =  T2  -  Tj.  The  slope  of 
Av/v  versus  surface  stress,  known  as  the  acoustoelastic  constant,  is  then 
obtained  at  various  frequencies. 

EXPERIMENTAL  RESULTS 

Figure  4  shows  Av/v  versus  surface  stress  at  3  and  5  MHz.  Figure  5  shows 
the  experimentally  obtained  acoustoelastic  constants  for  Rayleigh  waves  at  fre¬ 
quencies  between  3  and  5  MHz.  A  linear  least  squares  fit  of  the  acoustoelastic 
constants  to  Eq.  (12)  was  then  performed  to  obtain  experimental  values  of  the 
coefficients  80  and  Bj.  We  found  that  B0  =  0.009  1/GPa  and  =  -0.032  MHz/GPa . 

A  number  of  explanations  for  the  differences  between  the  theoretical  and 
experimental  determinations  of  B0  and  Bj  are  possible.  One  possible  explanation 
is  that  the  elastic  constants  used  in  the  calculation  are  not  representative  of 
this  sample.  This  specimen  was  cut  from  a  cylinder  that  had  previously 
undergone  an  autofrettage  process  (plastic  deformation  to  induce  compressive 
residual  stress  at  the  inner  surface).  Other  treatments  that  the  material  was 
subjected  to,  such  as  cold  rolling  and  forging,  can  also  have  drastic  effects  on 
the  elastic  properties  of  the  material. 

CONCLUSION 

As  expected,  the  relative  change  in  Rayleigh  wave  velocity  was  observed  to 
be  proportional  to  the  applied  stress,  with  a  proportionality  constant  known  as 
the  acoustoelastic  constant,  that  varied  with  frequency.  Furthermore,  the 
variation  of  the  acoustoelastic  constant  with  frequency  agreed  crudely  with  the 
Bq  ♦  B\/f  form  predicted  by  theory.  The  differences  in  magnitudes  of  B0  and  B^ 
can  possibly  be  attributed  to  changes  in  the  elastic  constants  caused  by  treat¬ 
ments  that  the  specimen  had  previously  undergone. 
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Figure  1.  Sample  and  transducer  arrangement. 
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SMCAR-AET-O,  BLDG.  351N  1 

SMCAR-CC  1 

SMCAR-CCP-A  1 

SMCAR-FSA  1 

SMCAR-FSM-E  1 

SMCAR-FSS-D,  BLDG.  94  1 


SMCAR-IMI-I  (STINFO)  BLDG.  59  2 

PICATINNY  ARSENAL,  NJ  07806-5000 


NO.  OF 
COPIES 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  1 

ROCK  ISLAND,  IL  61299-5000 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND,  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIB)  1 

WARREN,  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMAND 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 


DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLCBR-DD-T,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN :  AMXSY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5071 

COMMANDER 
HQ,  AMCCOM 

ATTN:  AMSMC-IMP-L  1 

ROCK  ISLAND,  IL  61299-6000 


COMMANDER 

US  ARMY  FGN  SCIENCE  AND  TECH  CTR 
ATTN :  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT'D) 


NO.  OF 
COPIES 

COMMANDER 

US  ARMY  LABCOM,  ISA 

ATTN:  SLCIS- IM-TL  1 

2800  POWDER  MILL  ROAD 
ADELPHI,  MD  20783-1145 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 

ATTN:  CHIEF,  IPO  1 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 
DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  MATERIALS  SCI  &  TECH  DIVISION  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  D.C.  20375 


NO.  OF 
COPIES 


COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 

ATTN:  AFATL/MN 

EGLIN  AFB,  FL  32542-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORA  ORY 

ATTN:  AFATL/MNF 

EGLIN  AFB,  FL  32542-5434 

METALS  AND  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  DIVISION 
505  KING  AVENUE 
COLUMBUS,  OH  43201-2693 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


I 


E: 


